Comment on "Quenches in quantum many-body systems: One-dimensional 

Bose-Hubbard model reexamined" 



O 

(N 

Oh: 

CD 
00 



o 

CD 

e 

I 



O 

in 

O 
0^ 
O 



X 



Marcos Rigol 

Department of Physics, Georgetown University, Washington, DC 20057, USA 
(Dated: September 3, 2010) 

In a recent paper Roux [Phys. Rev. A 79, 021608(R) (2009)] argued that thermalization in a Bose- 
Hubbard system, after a quantum quench, follows from the approximate Boltzmann distribution of 
the overlap between the initial state and the eigenstates of the final Hamiltonian. We show here 
that the distribution of the overlaps is in general not related to the canonical (or microcanonical) 
distribution and, hence, it cannot explain why thermalization occurs in quantum systems. 
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In Rcf. [l[, Roux argued that thermalization in quan- 
tum systems, after a quench, follows from the approx- 
imate Boltzmann distribution of the overlap between 
the initial state and the eigenstates of the final Hamil- 
tonian. Roux studied quantum quenches in the Bose- 
Hubbard model and concluded that, for small quenches, 
Cctl^ ~ K^'al^mi)!^ exhibits an exponential decay typ- 
ical of a canonical ensemble. The distribution of |CqP 
(or pm as denoted in enters into the computation 
the infinite-time average O of any observable O, where 
O = EalCa|'0„a, Oca = a.\d\^ . are the 

eigenstates of the final Hamiltonian, and \'ipini) is the 
initial state. This would explain why thermalization was 
observed in Ref. for small quenches. Here, we show 
that the distribution of |CqP is not related to the canoni- 
cal (or microcanonical) distribution and, hence, it cannot 
explain why thermalization occurs in quantum systems. 

We study a nonintegrable model of hardcore bosons 
(HCBs) in a linear chain with nearest-neighbor hopping 
t and interaction V , and next-nearest-neighbor hopping 
t' and interaction V' . We perform a quantum quench 
from two different initial states that have the same en- 
ergy Eq = {'^ini\Hfin\^ini), and hence the same effec- 
tive temperature T in the final Hamiltonian Hfin- 
We utilize full diagonalization to study eight HCBs in a 
24-site lattice with periodic boundary conditions. Trans- 
lational symmetry is used and the initial state is selected 
from the eigenstates of the initial Hamiltonian with total 
fc = 0. Further details about the equilibrium proper- 
ties and nonequilibrium dynamics of closely related HCB 
systems can be found in Ref. [31 and for spinless fermion 
systems in Ref. @. 

In Figs. [Ija) and [ijb), we show the distribution of 
\Ca\^ for the two initial states selected for our quenches. 
They are compared with the canonical weights corre- 
sponding to the effective temperature T. Two salient 
features can be seen in those figures, (i) The distribu- 
tion of ICqP exhibits a maximum around the energy 
_Eo = —4.62, while in the canonical distribution the 
ground state is always the state with maximal weight, 
(ii) The exponent of the exponential decay of |Ca|^ is 
different for the two initial states, i.e., it depends on the 
initial conditions and cannot be predicted by the effec- 
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FIG. 1: (Color online) Results for a quantum quench from: 
tini = 0.5, Vini = 2.0 (State 1), and tini = 2.0, Vini = 0.5 
(State 2), to tfin = 1.0 (sets the energy scale), Vjin = 1.0. 
In both cases t',^, = f'^,„ = 0.32, V/„i = V'f,^ = 0.32, Eo = 
—4.62, and T = 3.0. (a),(b) |CaP (points) and canonical 
(straight Une) distributions, (c) Energy distributions (d) n{k) 
of the two initial states, (e) n{k) of the infinite-time average 
for both initial states and n{k) of the standard statistical 
ensembles. Inset in (d), n(k = Q)aa vs E for all the eigenstates 

of Hfin- 



tive temperature of the system. From these results we 
conclude that, in general, the distribution of |Cap is not 
related to standard statistical ensembles and cannot ex- 
plain thermalization. 

The energy distributions p{E), where p{E) = (proba- 
bility distribution) x (density of states), corresponding 
to both initial states, as well as the one corresponding 
to the canonical ensemble, are shown in Fig. [He). They 
provide guidance to identify from which region of the 
many-body spectrum are the eigenstates that contribute 
to O, given the two initial states, and to the canonical en- 
semble result. As expected, all the energy distributions 
peak around Eq = —4.62 but they are all different from 
each other. We note that, as discussed in Ref. Q, the 
width of the energy distributions for both initial states 
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is expected to vanish in the thermodynamic hmit, as the 
width of the canonical distribution does. 

Figure [Iljd) depicts the momentum distribution func- 
tion [n(fc)] of both initial states, which are clearly dif- 
ferent from each other. In Fig. [ije), we present the 
results of the infinite-time average of n{k) for the two 
initial states, and compare them with the predictions of 
the microcanonical and canonical ensembles. The time 
averages are virtually indistinguishable from each other 
(independence of the initial conditions) and agree with 
the microcanonical prediction, i.e., thermalization takes 
place. The canonicalprediction is slightly different due 
to finite size effects [J, Q . Considering that | Cq | ^ ver- 
sus E is different for both initial states and from the 
microcanonical weights, one can understand thermaliza- 
tion in terms of the eigenstate thermalization hypothesis 



(ETH) 0] . ETH states that the eigenstate expectation 
value of generic few-body obscrvables are very similar be- 
tween eigcnstates that are close in energy. From ETH it 
follows that no matter the weights one uses in the av- 
erage over a narrow window of energies, the result will 
always be the same. The validity of ETH for n{k)aa 
around Eq = —4.62 is corroborated by the behavior of 
n{k = Q)aa versus E in the inset in Fig.[Hd). 
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